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Calculation of Thermodynamic Functions for 
Polyatomic Molecules' 

Harold W. WooJley 

Formulas are given from which the thermodynamic functions may be obtained for 
polyatomic molecules not exhibiting special phenomena, such as internal rotation. The 
effects of first and second order anharmonicities and rotational-vibrational interaction are 
treated in detail, as also are the angular momentum effects of doubly degenerate vibrations 
in linear molecules. This paper presents formulas for the corrections carried to higher order 
than has hitherto been available. 

1. Introduction 

There is occasionally need for accurate thermodynamic functions for polyatomic molecules. 
A definite calculating procedure exists for approximate calculations via the rigid-rotator har- 
monic oscillator treatment, with small corrections for effects of anharmonicity, stretching, etc. 

This paper presents formulas for the corrections carried to higher order than has hitherto 
been available. Calculations can thus be made via the approximation formulas with accuracy 
probably comparable with that of the direct summation method using extrapolated formulas 
for the energy levels, aside from effects related to vibrational and rotational cut-off ; which have 
been omitted. 

2. Nonlinear Polyatomic Molecules 
The partition function for the internal molecular energy states may be written as 

Q=QbQrcQh.O.Qc (i) 

For nonlinear molecules, 



Q R =(lh)^(TlA B C )(kT/hcy (2) 

is the rigid rotator partition function for the lowest vibrational energy state, where A 0y B , and 
Co are the rotational constants for the ground vibrational state, and y is the symmetry number. 
Qh.o. is the partition function for vibrational levels of harmonic oscillators. Q Rc is a correction 
factor including (1) a factor 1 + pT. . . for rotational stretching, as shown by Wilson [l] 2 and 
(2) a factor l+0 1 /T+0 2 /T 2 for a low temperature quantum correction for rotation as shown by 
Stripp and Kirkwood [2]. 3 Q c is a factor including rotational-vibrational interaction and the 
anharmonicity corrections, now to be treated in detail. 

The rotational-vibrational interaction in its simple form arises from the vibrational de- 
pendence of the analogue of eq (2) for higher vibrational states. Before summation over 
vibrational quantum numbers, the dependence is introduced by including the factor 

#t>=exp[2a^+2a^+2a<^ (3) 

equivalent to a form introduced by Benedict [3]. 4 With his formulation, the constants in the 
factor R v come from sums of corresponding constants for A VJ B v , and C v . A summation pro- 
cedure adequate to treat this with the anharmonicity effects follows, extending the latter 
beyond the linear terms of Stockmayer, Kavanagh, and Mickley [4]. 

i Related to a thesis submitted to the Graduate School of the University of Michigan in partial fulfillment of the requirements for the Ph. I ). 
degree in Physics. 

2 Figures in brackets indicate the literature references at the end of this paper. 

3 Expressions representing the T 2 rotational stretching effect and the T -3 low temperature correction as well as a low temperature correction 
for symmetry effects for different classes of rotat ional levels for the asymmetric rotator are included in the thesis. 

4 An equivalent representation in slightly different form is given in the thesis. 
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For a given vibrational level, specified by vibrational quantum numbers v i} l<i<n, the 
statistical weight is the product of the n quantities 

Pi=(v i +d i -l)\/[v i \(d i -l)\] ) (4) 

where d t is the degeneracy associated with the ith. fundamental frequency. 
The summation over the vibrational levels may be indicated as 

Q=QrQ*c?Rv ex V [-hcG v /kT]U Pi (5) 

with the vibrational energy in wave numbers given by 

G 9 =2o} i v i +2x ii v i (v i — l)+Zx ij v i v j +^y iii v i (v i — l)(v i —2)+2y iij v i (v i —l)v j +2y m v i v j v k . (6) 

The relations between these constants and the more usual ones are indicated in appendix 1. 

The Boltzmann factor is separated into two factors, one for harmonic oscillators fitting 
the lowest levels and a second giving the alteration due to anharmonicity 

exp[-Ac67^T]-exp[-k2co^^T]exp[-Ac(^-2co z ^)/^r]. (7) 

Q vc , the product of the second factor on the right hand side of (7) with R v , is expanded into 
a power series whose terms are polynomials in the v'a. Defining a function 

f t (Z)=v t \l(v t -Z)\ (8) 

it can be shown that 5 

2> <e -« i »,:=(l-e-« i )-<S (9) 

where u t =hcw t /kT, and that 

J2f t (Z)p l e-^i = [(d t +Z-l)\I(d i -l)\]-[e- Zu iKl-e^iYi+ z ]. (10) 

Vi = 

The factor Q vc can be expressed as a sum of terms, each linear in whatever //s it contains. 6 
The results of multiple summation over all #/s follow by application of eq (10). If the right- 
hand side of eq (10) be designated by /*(Z), one has 2 v Q v Q vc =Q Ht0 ,Q c , where Q v =exp [—2 i u i v t ] 
Qh.o. = ^vQv and where Q c is the same function of f*(Z) as Q vc was of /*(Z). Thus eq (1) 
has been evaluated, Q H .o. being given by 

&.o.==n(i--«-«o-*. (ii) 

i 

There is some advantage in using the logarithm of the partition function from this point 
on because of analytic cancellations of some higher order terms. The terms removed are cross 
products without linkage in the subscripts, such as x tj x kh x u x jk , and a 2 % 2 , for example. 

The contribution to In Q c due to x H and x tj terms, with r t =e~ u i and s i =(l — e~ u i)~ 1 is 

+a£i<Wjr«r,[l+d«r,+d^]^ (12) 

with terms of higher order for this and other contributions listed in appendix 4. At high 
temperature the first term varies as T and the second as T 2 . The terms in the appendix are 
higher than T 2 in principal dependence at high temperature. 

5 See appendix 2. 

e Relations giving the simpler polynomials in terms of the/'s will be found in appendix 3. 
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Similarly, the contribution of higher order anharmonic terms is 

- (hc/kT)J2 { yu4i(di+ 1 ) (d i +2)r*s!+y Uj d i (d i +l)d j r 2 i r J s*s j +y m d i d j d k r i r j r^ } . (13) 

The contribution due to rotational -vibrational interaction is 

y2 { a i diriS i +a ii d i r i (l+d i r i )s*+a ij d i d j r i r j s i s j J r ^oid-r^ } (1 4) 

as far as a's alone are concerned, while their mixing with anharmonicities gives 

— (hc/kT)^2 { 2a i x ii d i (d i + 1 )r?s? +a i x i jd4jr i r^8 j } . (15) 

In all of these (multiple) summations, if the terms are completely symmetric in the indices 
that are present, such as i, j, etc., the symmetric indices are to progress in size as with i<ij<ik, 
and so forth. For terms not completely symmetric in the indices, the indices not symmetric 
are simply not to be equal, thus i^j^k, and so forth. 7 

3. Linear Molecules 

There is an additional kind of vibrational-rotational interaction, concerning vibrational 
angular momentum, which is important for linear molecules. When doubly degenerate vibra- 
tions occur, the vibrational energy, eq (6), contains an additional term g=2 g^lilj. For 

i<j 

linear molecules, in addition, there is dependence of the rotational energy upon / where l=2l { , 
giving for the two together 

F{J, l)=g+B v [J(J+l)-l>]-D v [J(J+D-l 2 ] 2 , (16) 

where J— I, Z+l, Z+2 . . ., with double degeneracy for |Z|^>0. 

The sum over the rotational levels for J> \l\ follows from a half integer summation formula 
of the Euler-Maclaurin type [5] 



n=ri\ J n 



giving 



and 



/(^7i+(l/24)r(^ 1 -|)-(7/5760)r / '(n 1 -|) + (31/967680)/^(7i 1 -|) . . .,(17) 



S (2J+l)exphcB[J(J+l)-l 2 ]/kT=kT/hcB+l/Z + . . . (18) 

J=\i\ 



Q= S (2J+l)e- hcF/kT =(lly)[kT/hcB+l/S+2(kT/hcB) 2 D/BYe- hcg/kT . (19) 

The exponential in g is expanded as a power series and the summations made with v t increasing 
by steps of two from l t to oo f then summing l t from — oo to + oo . If g is extended to include 
possible contributions to In Q c varying as T 2 at high temperature, while retaining the g tj term 
needed for nonlinear molecules, as in 

g = ^2gJi + S#*M + iLjQuRvi + ^0,gukliV k , (20) 

i KJ 

the contribution of g to In Q c becomes 

-2(/^T)2{^^+^^ 

+ r?)s?+20? y r<r^s?+40„3« < r^ (21) 



7 For a case with 3 frequencies, for example, xuxa which is not completely symmetric in the indices would include x\\X\z and also X33X31. But 
xaXjjXij has complete symmetry in the indices, and £33X11X13 is sufficient without £33X11X31 also. 
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Constants of the latter two types in eq (20) are not experimentally available at present. For a 
linear molecule, when there is only one doubly degenerate frequency, with l=l 2 , the —B v l 2 
term of eq (16) is commonly absorbed into the g^% term. In that case, g H of eq (21) repre- 
sents g22~\-B v of the spectroscopist. 

The free energy function for a mole of a substance in the standard ideal gas state is given by 

-(F o -E° )/RT=\n Q+(5/2) In T+(3/2)ln M+ln (2tt /iV) 3/2 & 5/2 A 3 ^o, 

P being atmospheric pressure, M the molecular weight, and the others known physical con- 
stants. The other functions follow by differentiation as given in appendix 5. 

A tabulation by Pennington and Kobe [6] may also be used for convenient evaluation of 
the more important of the contribution here treated. 
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5. Appendices 

5.1. Appendix 1 
The relation of the constants of eq (6) with those in 

0.=S««(»«+d«/2)+:E^^^ 

+ J2yw(v i +d i /2nv j +dj/2)+ Ti yij*(Vi+di/2)(vj+d } l2)(v k +d k /2) 

i^j i<j<k 

is 

^=^,+ (^+l)^+SH^o-+(l/4)(3^ + 6^+4)^ 

+S(i)(*+W«(+S(W(»+SS(W«i 

3 j j<k 

Xu=x ii +(i)(Sd i +d)y iii +^,(d j /2)y iij 

j 

x ij =x ij + (d t + l)yuj+ (dj+ 1 )y w + 2 (d k /2)y ijk 

k 

yiii—ym, yuj—yaj & n d yijk—yw- 

5.2. Appendix 2 
To obtain eq (9) and (10), begin with the identity 

Multiply by x z ~ l and differentiate n times with respect to x, getting 
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For the second member of the identity, let x—l=t, giving 

W®+W or ^[g^R^- 1 ] 

(I 71 

For ">Z— 1, the only contribution is from the one term with r=Z, giving -j— t~ l or (— l) n n\ 

x l-Z+n 

(x— \)~ l ~ n . Multiplying both members of the identity by (— l) n , _ 7 w > one has 

(v+n-Z)\x- v _ n\ z _ r -iw-« 

£i{v-Z)\(n-Z)\ (n-Z)\ x U x } ' 

Inserting e~ u * for x, one may obtain eq (10) by using n=Z J rd i —l, with eq (9) as the special 
case given by Z=0. 

5.3. Appendix 3 
Functions of v in terms of (f(s)=v\/(v—s)l: 

« 2 =/(2)+/(l) 

D»=/(3)+a/(2)+/(i) 

* 4 =/(4)+6/(3)+7/(2)+/(l) 

» 6 =/(5) + 10/(4)+25/(3) + 15/(2)+/(l) 
fl 6 =/(6) + 15/(5)+65/(4)+90/(3)+31/(2)+/(l) 

e»-t>»=/(3)+#(2) 
B»-t^=/(4) + 5/(3) +4/(2) 

» 4 -2« 3 + « 2 =/(4) +4/(3) +2/(2) 
fl 5 -2y 4 +«; 3 =/(5) + 8/(4) + 14/(3) +4/(2). 

5.4. Appendix 4 
Additional terms in the formulas given earlier are as follows: 
Equation (12): 

-(l/6)(/W^r) 3 2:{4a; 3 i J i (rf i +l)rHl + (8(/ i +12)r t +(8^+22d i +15M+(2rtT+4rf i +2)rfK 
+x?,<M/i»v[l + (3<*i+ 1 )r t + (M j +l)r j +dM+dW+(Qd l d i +3d t +3d j +l)r i r } +d 2 t ijrj 

+ (2d<+ l)r i +2d j r j +djr i r j }sts 2 + 3x 2 i x ik d i d j d k r i r j r k [l + (2di+l)r i +^+d/ < ry]sfs?8» 
+24a; 1J ^/r 4 /? 1 (d < +l)^(dj+l)rfrfs?sf+12x«a; < /i; ( t(i < (d < +l)d^j ; rfr/it(2+r J )sJs^ J | : 

+ \2x H x tj x jk d i {d i +\)d 1 d t r^rf k s 3 i s 2 s k + &x ii x ik x tl d l dfl k d l r i r j r k ri(l +r i )sfs i 8 i 8, 
+ Qx tj x {k x j id i djd k dir i rf k riS 2 s 2 s k Si + &x iJ x tk x jk d i d ] d k r i r j r k (\+dir i +d j rj+d k r k )s 2 s 2 sl} 
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Equation (13): 

-(fc/i:nS{W<(d<+i)(*^ 

+ z im d i (d i + 1 )d j d k ?> 2 r j r k s 2 s j s k + z iijj d i (d i + 1 )d j {d j -\- 1 Y i r*s*s 2 j + z im d i d j d k d l r i r j r k r l s i s } s k s l } 

+ (l/2)(hclkT) 2 j:{12x u y u Md i +l)(d i +2)rni + (d i +iy 

+2(d i +l)r i ]sts j +4:X u y jji d i (d i +l)d j (d j +l)r 2 ^ 

+ 6x ij y iii d i (d i +l)(d i +2)d j rlr j sis j +2x ij y iij d i (d i +l)dffi 

+4:X ij y U kd i (d i +l)d j d k r 2 r j r k s!s j s k +2x i3 y kki d4 J d k (d k ^ 

s j s k Si+2x ij y m d i d j d k r l r j r k (l+d j r i +d J r j )s 2 s 2 s k } 

Equation (14): 

+dffl]$$+a ii $4jr i r j (l+d i r i )s?s j +a ijJ 4id j d k r i r j r^ 
Equation (15): 

-(Ac/& 7)2(0^^(^+1)^ 

+ (2d*+ 1 )rjsf^+a^^^r^[l +d i r i +d j r j \s 2 s 2 +a ij x kj d i djd k r i r j r k s i s 2 s k } 

+ UhclkT) 2 2{4a i x 2 J i (d i +l)r 2 [l + m^ 

+ (2d i +l)r i +d j r j +d j r i r j \s*s 2 +4a i x ii x ij d i (d i +l)dffi 

+ l)d j r 2 i r j s*s 2 +2a i x ij x ik d i d j d k r i r J r k (l+r^ 

Equation (21): 

-(l/3)(hc/kTyjl{gUr i +50r 2 i + lMrl + 
+3x ij g 2 li d j r i r j (l + 19r t + lQr 9 i +fi)^s,+12x iJ g 1 ig i fi'^ s (l +r t ) (1 +r i )* l «»J 
+123^^(4+55^+7071+9^ 

+3a^«<Z i 7'<r>[l +8r i +5r?+d J (ry+4r 4 r y +r 2 < r i )]s 4 ^+ 12» w a: i ^ i< d j (d i + l)r^J(l +*•<)*«*' 
+Sx ij x ik g ii d j d k r i r j r k (l +4:r i +r 2 i )s\s j s k +6x ij x jk g i 43dkrir j r k (l+r i )ey j s k } 

5.5. Appendix 5 
For a term of In Q of the form 

z=vu?f 1 llr? t 8? i =Wf u 

i 

with /i = l + S'i r ? i an d u=hc/kT and with #,. p, m { , %, t it and g* constants, the contributions 

to the various' thermodynamic functions are: For -(F°-E° )IRT, z^Wfc; for (H°-E°)IRT, 
—udz/du=Wh 2 ; for C° p /R, u 2 d 2 z/du 2 =Wh 3 ; with h 2 =figi+f 2 , h=fi92+2f 2 gi+f 3 , where 

i 

^ 2 =— p + 5>^? r A ? +0i. 

i 

Washington, January 11, 1956. 
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